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Abstract 

The nucleus-nucleus Glauber amplitude in the tree approximation is studied for 
heavy participant nuclei. It is shown that, contrary to previous published results, 
it is not unitary for realistic values of nucleon-nucleon cross-sections. 

1 Introduction 

Scattering on the nuclei is commonly studied in the Glauber approximation, which can 
be rigorously derived in Quantum Mechanics provided the transverse momenta trans- 
ferred to the projectile are much smaller than its longitudinal momentum. With certain 
reservations it can be generalized to the high energy region where the elementary nucleon- 
nucleon (NN) amplitudes become predominantly inelastic. For the nucleon-nucleus (NA) 
scattering the Glauber approximation has a transparent probabilistic interpretation. If 
the target nucleus is heavy, with atomic number A » 1, the Glauber formula acquires a 
simple eikonal form, which clearly shows that the resulting amplitude is unitary, that is 
its modulus is smaller than unity at fixed impact parameter. 

With the advent of collider experiments nucleus-nucleus (AB) scattering becomes an 
important physical object. The Glauber approximation can be easily generalized to the 
AB case and it was in fact done very long ago. The Glauber formula for AB scattering 
looks very similar to the hA case. At fixed impact parameter b the scattering matrix S 
is assumed to be a product of nucleon-nucleon scattering matrices s averaged over the 
transverse distributions of nucleons in both nuclei: 

S(b) = (f[f[s(b-x l + x' k )) AB , (1) 

i=l k=l 

where Xj and x' k are the transverse coordinates of the nucleons in the projectile and target 
nuclei respectively and in absence of correlations in both nuclei averaging < ... >a,b 
means 

- A B 

(F(x i ,x' k )) AB = J ]J]Jd%d 2 x' k (T A (x l )T B (x' k -b)F(x 1 ,...x A ,x' 1 ,...x' B )). (2) 

i k 

Here T A (x) and Tb{x') are the standard nuclear profile functions normalized to unity. 
However in contrast to the NA case the content of the Glauber formula for AB scattering 



a 



b 



c 



Figure 1: Examples of disconnected (a), tree (b) and loop (c) diagrams for the Glauber 
AB amplitude with A = B = 2 



turns out to be much more complicated. Presenting in the standard manner the NN 
scattering matrix 

s(b) = l + ia(b), 

where a is the NN scattering amplitude, one obtains from a set of terms corresponding 
to different ways the nucleons from the projectile and target may interact with each other. 
Each of these terms may be illustrated by simple diagrams indicating these interaction. 
Some examples are shown in Fig. 1. for two pairs of interacting nucleons in the pro- 
jectile and target. One observes that in contrast to hA case the diagrams may contain 
disconnected parts (FigJTJa) and, most important, loops (Fig{TJc), which involve internal 
integrations over transferred transverse momenta and thus NN amplitudes for non-zero 
transferred momenta. 

Loop contributions depend not only on the total NN cross- sect ions as the tree diagrams 
but also on the form of the differential NN cross-section. Their calculation is difficult and 
unreliable, so that in most applications to heavy nuclei (A, B » 1) loop contributions are 
simply neglected. The typically used approximation is the so-called optical approximation, 
which corresponds to taking into account only the simplest contribution (a single NN 
interaction) for each connected part. In the optical approximation 

S°?\b) = e F °^ b \ (3) 

where the optical eikonal F opt (b) is 

F opt {b) = iaAB J d 2 xT A {x)T B {b - x). (4) 

Its advantage is simplicity and obvious unitarity. The natural question, which has been 
long discussed in literature, is the quality of the approximation which neglects loops (the 
tree approximation) for the case A,B » 1 and in particular its relation to the optical 
approximation. In PQ in the limit A, B » 1 a closed formula was obtained for the tree 
approximation to the Glauber amplitude, different from the optical approximation but 
also unitary in the above mentioned sense, that is with a modulus smaller than unity at 
fixed impact parameter. Their result gives the tree eikonal in the form 



F®(b) = - [ d 2 xf( 7A , lB ), 
a J 



where 

ja(x) = —iaAT^ix) , ^b{x) = —iaBTs(b — x), 

k 

/(7a, 7b) = + v i + u ^) - 7a - 7s 

i=i 

and and t>/ are the k solutions of the transcendental equations 

v = 7 B e"". 

The number of solutions may be k = 1 or k = 3. In the latter case it is assumed that 
U\ > U2 > U3 and v\ < v% < V3. One can check by numerical studies that the resulting 
F(b) is always negative, so that \S(b)\ < 1. However later in [2] it was claimed that in fact 
at A, B » 1 the sum of all tree diagrams is exactly given by the optical approximation. 

The aim of the present study is to resolve this contradiction for the asymptotic of 
the tree approximation to the Glauber AB amplitude for large A and B. Our results are 
different from both p] and [2]. Unfortunately they are also much gloomier. Namely 
we find that at A, B » 1 the sum of all tree diagrams inevitably becomes non-unitary, 
so that taking loops into account is absolutely necessary for the physically meaningful 
amplitude. Thus the optimistic hope that the dependence of the AB amplitude on the 
behaviour of the NN amplitude at non-zero momentum transfers is hardly probable. [1] 
is unfortunately not justified. The reason why our results are different from the previous 
ones lies in the details in which simplifications related to the asymptotic at A, B » 1 
are made, in particular in the not completely accurate use of the saddle point method in 
the previous derivations. As we shall see apart from the saddle points which lead to the 
quasi-optical approximation there exist others which may give the dominant contribution 
and destroy unitarity. 

Our study will be based on the equation which expresses in a compact manner the 
sum of all tree diagrams, obtained by one of the authors some 30 years ago [3] . The good 
quality of this equation is that it is valid for arbitrary finite (and even small) values of 
A and B and therefore presents an adequate starting point to investigate the asymptotic 
at large A and B. To simplify we shall limit ourselves with the case of collision of two 
identical nuclei A = B. Since the problem we address does not depend on the form of the 
transverse distribution T^{x) we further simplify our study by assuming that T(x) does 
not depend on x inside the nucleus: 

Ta(x) = 9{Ra — \ X \)~T>2 (5) 

where Ra = A l ^R Q is the radius of the nucleus. 

2 The tree amplitude for AB scattering in the Glauber 
approximation 

2.1 General A, B and b 

In [3] an expression was derived for the sum of all tree diagrams for the S matrix in 
the Glauber approximation to nucleus-nucleus scattering valid at arbitrary finite atomic 
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numbers A and B of colliding nuclei. At a given impact parameter b 

S(h) - AWl I dTdT ' J(T+T')-z ib ,ry) (6) 

Here 

Z = -ij (d 2 x)w(irT A (x),ir , T B (x')) (7) 

and 

(d 2 x) = d 2 xd 2 x'5 2 (b — x + x). 

As indicated, a is the nucleon-nucleon forward scattering amplitude and Xa(x) and Tg(x') 
are the nuclear profile functions at transverse coordinates x and x' respectively. W is the 
effective classical action for the effective quantum theory of two fields u and v 

W(p,p') =y(u(p,p')Xp,p')p,p',), 

where 

Y = -vu - ip(e~ u - 1) - ip'(e' v - 1) 
a 

and u and v satisfy a pair of transcendental equations 

u = —iap'e~ v , v = —iape~ v . (8) 

Here p = itTa(x) and p' = zt'Tb(x'). 

A detailed derivation of this formula can be found in [3]. For convenience it is briefly 
reproduced in Appendix. In this section for the simplified case of constant T4 and Tg 
inside the nucleus we transform this formula to the form suitable for our analysis at large 
A and B. 

It is trivial to see that W is different from zero only in the overlap region. Indeed for 
Ta = and so p = we have v = and for Tb = and so p' = we have u = 0. In both 
cases Y = 0. Therefore integration over x and x' is extended over the overlap region only. 
For constant Ta,b it gives precisely the area of the overlap region G(b). Taking this into 
account and expressing W via Y we find 

K ) ~ A-K 2 l A+B + 2 J T A+l T lB+l ' 



where 
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and equations (jHJ) become 

u = K'T'e' v , v = KTe~ u . (11) 
Separating in the exponent in ([9]) the terms proportional to r or r' we rewrite Eq. ([9]) as 



To avoid solving transcendenal equations ffTTT) we pass in (fT2j) to the integration over k 
and v , since it is trivial to express r and r' via w and f from ffTTT) but not vice versa. To do 
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this we have to find the Jacobian. We denote £ = kt and 77 = k't' Direct differentiation 
of ( jTTj) gives 

du dv du u dv dv v du dv du 
<9£ c>£' dr] r] dr] } <9£ £ <9£' <9?7 9ry 

From these equation we immediately obtain 

du u dv v du uv dv uv 



dr] 77(1 — uv)' <9£ £(1 — uv) 7 <9£ £(1 — -ut>)' <9?y 77(1 — uv) 
As a result we find the Jacobian 

d(u,v) uv 



J 



d(tv) £77(1 -uv)" 



At small x and ?/ we obviously have u ~ 77 and t> ~ £. So choosing the initial contours 
in x and y around the origin small enough we find that integrations over u and v will go 
also around small contours around the origin, which can then be transformed unless we 
come across some singularities in u and v. Expressing r and r' in terms of u and v as 



1 „ / 1 t 
r = -i>e , r = —ue 

K K 



we transform Eq. ( fl2l) into 



This formula is the starting point for our investigation. 
2.2 Case A = B and 6 = 

Our formula for S(b) greatly simplifies in the case of central collisions of identical nu- 
clei, when A = B and 6 = 0. In this case the complicated exponents in the first two 
exponentials in ffTBj) are absent and we find a simple expression 

m = i^SBS""/ (|pr<l - »)e^ e -*<-^'. (14) 

It is straightforward to find this S matrix in the form of a finite sum of terms. Integrations 
over u and v obviously give the coefficient before term (uv ) A in the expansion of the rest 
part of the integrand in powers of u and v (with factor (27ci) 2 which cancels with the 
analogous factor in front of the whole expression in (1141) ). So our problem reduces to the 
expansion in powers of u and v of the three exponentials in (fl4"l) . We find 
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e \ in -y- e i»w 



1 / l\ n3 /l \ 111+112 

V ———-( — ] [ A) u ni+m v n 2 +n 3 



ni i« 3 ri x \n 2 \n z \ \ ikJ \ in 
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t — v uniyn2 ( 1 A \ ni+n2 ~ 2n (15) 

Integration over u and v gives the term with n\ = n-i = A. Without the factor (27ri) 2 and 
the one in front of the whole expression (fT4l it is 

a 1 / i \"/ i \ 2(A ' n) 

^ n!p-n)!] 2 V W V m ~ J 

= E UM m ( - 1 ) ( J - ^ 2(A "" } - ( 16 ) 

^ n\[{A-n)\] 2 V m7 

7 = -iAk. (17) 



where we defined 



From this expression one has to subtract the second one which comes from the term —uv 
in the Jacobian. Obviously it gives the (A — l)th term in the expansion of the integrand 
in powers of u and v and is obtained from ffl6|) by the substitution A — > A — 1. 
Collecting all the factors we get 



S(0) = (A!) 2 (l- 7 )^{£ 



1 ( 7 



n=0 



n\[(A-n)\} 2 \A(1 - 7) 1 



^ 2 (1 - 7) 2 ^0 " 1 " n W \ A ( l ~ 7) 2 , 

Note that in fact the scattering amplitude a is pure imaginary at high energies: 

i 

a = -a, 

where a is the total cross-section for pp collisions. Threrefore 7 is positive and so each 
term in the two sums in ( fl8l) is positive. This expression gives a simple closed form for 
the 5* matrix for collisions of identical nuclei at b = 0. The term with n = in the first 
sum is independent of the scattering amplitude a and equal to unity. The rest terms give 
the amplitude with factor i. 

Note that (1181) contains powers (1 — 7)™ = (1 + Ain) n . This is the origin of difficulties 
related to AB scattering. A similar formula for NA scattering contains powers (1 + in) 71 . 
Since n is small, of the order A~ 2 ^ 3 , summation over n does not violate unitarity. In 
contrast, for the nucleus case n is substituted by An, which grows with A as A 1 / 3 . As a 
result factors (1 — 7 ) n grow like A n l 3 at high A and n and as we shall discover make \S\ 
also grow. Unitarity remains valid only for values of 7 just slightly above unity, which is 
certainly not satisfied at high enough (and physically interesting) A. 

Eq. (118p makes it feasible to perform numerical calculations of the tree Glauber 
amplitude, since each of the two sums contains only positive terms. Numerical results 
demonstrate that S(0) is unitary, that is 15(0) | < 1, at any value of A only provided 
< 7 < 1.42. The last condition means that 



(19) 



To compare, a similar condition for hA scattering 

° <1 
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is always satisfied for large A. However in our case, with R A = RqA 1 / 3 ([in]) reduces to 



o-p-2 <L42 - ( 2 °) 
With cr ~ 27t_Rq it is always violated at physically relevant large A. 

3 The asymptotic of the amplitude at b = and A = 

B oo 

In this section we shall derive analytic asymptotic formulas which agree with our numerical 
results. 

We rewrite (fill) as 

where 

P(u, u) = — ln(-uu) — it — u H — (m> + u + u ). 

7 

We estimate the asymptotic by the saddle point method. The saddle point is determined 
by the equations 

P u = — - 1 + -{v + 1) = 0, 

u 7 

P v = l + -(u + l) = 0. (22) 

f 7 

The second derivatives are 

P =1 P =1 P = l - 

it f z 7 

Eqs. (1221 have two symmetric solutions 

u% = Vi = 7, m 2 = ^2 = -1- 
1. Saddle points w s = t> s = 7 

Passing to variables u — 7 = z£ and u — 7 = ^77 we rewrite the integral in (1211) in the 
vicinity of the saddle point as 



\2 „,2A 

. e AP(y,7) ( 

Subsequent actions depend on the sign of eigenvalues of the matrix 

m = — ( 1 Tl 



The eigenvalues are 

2^(1 + 7) and ^(1-t)- 

We have to consider two cases 7 < 1 and 7 > 1. 
1.1 Case 7 < 1 

In this case we can safely extend the integration regions in both £ and 77 to the whole 
real axis to find 

We further use 
and 

to finally find 



A\ = A A e~ A V2^A (24) 
AP(7, 7) = -A In 7 2 - A7 + 2 A 



S(0) = y]l- 7 2 e-^ 7 . (25) 

1.2 Case 7 > 1 

In this case we have to rotate one of the variables £ or fj which diagonalize matrix M 
which corresponds to eigenvalue 1 — 7 by angle ±7r/2. Then we get instead of ( 1231) 

1 = ±1- 



and for 5(0) 



S(0) = ^iyy - le~ A7 (26) 
2. Saddle points u s — v s — — 1 

In this case the prefactor (uv ) _1 — 1 vanishes at the saddle point and we have to study 
it in the vicinity of the saddle point. We put u = — 1 — z£ and v — — 1 — ir) to find 

1 = (1 - i£ - £ 2 )(1 - irj - rf) - 1 = -i£ - irj - £ 2 - r/ 2 - £77. 

Since the leading term vanishes we have to expand the exponent up to terms of the third 
order in £ and 77 

P(u, v) = P(-l, -1) - V + r? 2 ) - ifr + + r/ 3 ), 

2 70 

so that in the vicinity of the saddle points the integrand in fl2~T]) can be presented as 
_L _ \ \ e AP{u,v) _ _ e A(P(-i -V-i{?+^nhW)Q(t^\ 



-uv 

where the polynomial Q(^,rj) is 



So we find 
where 

J = JO ~ g-Ul +J2) 

with 

/ 9 1 d s . _ / d 2 1 «9 2 x . _ ,1 <9 2 X 

Jo = / didrte-^ 2 ^- 2 ^ 

and 

In the exponent the matrix in variables £, 77 is now 

Consider the case 7 > 1 Then a > (3 and both eigenvalues of matrix M are positive. 
Then we immediately get 

7r . _2tt7 2 (27 - 1) . _ 127T7 5 . 12vr7 4 

° ~ ^a 2 -(3 21 J ° ~ A 2 ( 7 2 - l) 3 / 2 ' Jl ~ A 3 ( 7 2 - l) 5 /2' j2 ~ ~^3( 7 2 _ 1)572 

Using also 

AP(-1,-1) = 2A- - 

7 

and the asymptotic f[2"4"j) we finally find for 7 > 1 



5(0) = -— — ^ TTT^e^ 21117 "^. (27) 

v 7 A(7 + 1) 5 / 2 (7 - 1)V2 v 7 

So for 7 > 1 5(0) is always negative and its modulus exponentially grows with A unless 
7 2 < exp(l/7), that is 7 < 70 = 1.4215 when it exponentially falls. 

The asymptotic for 7 < 1 requires additional rotation in variables which diagonalize 
matrix M. Up to its sign it can be found just by the analytic continuation from the case 
7 > 1. So for 7 < 1 

m - ^ (7+1 )4-7)'/^ <21 " 7 ""' (28) 

One observes that for 7 < 1 the leading contribution comes from the saddle point 
u s = v s = 7 and for 7 > 1 from the saddle point u s = v s = —1, the latter contribution 
restricting the region of 7 > 1 where unitarity is fulfilled 



4 Non-central collisions of identical nuclei 

We introduce 7 according to (JT7j) and put 

G{b) = X{b)7rR 2 A , < A < 1 (29) 

to rewrite Eq. ffT51) as 

*W = ^/^-lK<-> (30) 

where now 

P(u, t>) = — ln(ut> ) — m — v H — (w + it + w ) H (t>e" + ue'"). 

7 7 v 

The derivatives are 



Pu = - 
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it 7 
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Pv = ~ 


- -1 + - 






v 7 
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1 - A 



7 



(we" + e u 



P, tl , = — 2 + ^—ve u , P uv = - + ^—{e u + e v ).P vv = - + ^-^ue v , 
u 7 7 7 ^ / t> z 7 

As before we seek for symmetric stationary points u s = v s . Then we obtain an equation 

1 A 1 - A 
— + - + 

u s 7 7 

We have the same solution w 2 = v 2 = — 1 and a new solution «i = v\ which satisfies 

1 A 1 - A 



u s + l)( + - + e Us =0 



e Ul = 



Ml 7 7 

or 



Ml = A + (1-A)e^ (31) 

The actual value of u\ for given A and 7 can only be found numerically. 
1. Saddle points u s = v s = u±. 
We find at the stationary point 

P = -2 In mi + 2 - 2ui + -«?, 

7 

J_ _ A _2_ A 

-'tin Pvv 9 "I - 1 Wlj Puv • 

Uf 7 Mi 7 

The determinant of the quadratic form in w, t> in the exponent turns out to be 



1 1 



It is positive for u\ < 1 and arbitrary A < 1. So for u\ < 1 we find the asymptotic 
5(6) 



- J 4(2n 1 -A u 2 +21n 2L 



^ l_ U 2(l_ UlA/7)S 



(32) 



If b = and A = 1 then u% = 7 and this asymptotic passes into ( 1251) . The bracket in 
the exponent in (1321) is always positive and diminishes with A, which implies that the 
asymptotic gets less falling with the growth of b. At b = 2R& and A = 0, when the nuclei 
only touch each other, the exponent vanishes, which corresponds to S(b = 2Ra) = 1 as it 
should be. 

If Mi > 1 then the asymptotic can be obtained by analytic continuation of (l32il . It 
remain to be falling with A. However in some regions of A and U\ it becomes pure 
imaginary. 

2. Saddle points u s — v s — — 1. 

At this saddle point we find 

7 e 7 

Together with the prefactor in (1301) it gives an exponential factor in the asymptotic 

a(21h 7 - a_ aa) 

p \ T e T / 



(33) 



It infinitely grows at 7 > 70(6) where 



2 

2 7o In 70 = A + -(l-A . 

e 

The value of 70(6) steadily but slowly grows with A (see Table 1 ). These values restrict 
the region in which the amplitude remains unitary at b > 0. 



Table 1. 70 as a function of overlap A (b) 



A: 





0.2 


0.4 


0.6 


0.8 


1.0 


7o: 


1.3211 


1.3416 


1.3620 


1.3820 


1.4019 


1.4215 



Note that at b = 2Ra and thus A = the exponential factor (133|) grows with A unless 
7 < 7o(0) = 1.3211, in spite of the fact that S(b = 2R&) = 1. This seeming contradiction 
is resolved due to vanishing of the prefactor at exactly A = 0. As soon as A gets slightly 
greater than zero, the S'-matrix becomes very large and negative for 7 > 70(6). For 
example for A = 0.001 and 7 = 2.0 one finds S(b) = —0.174 ■ 10 5 , which illustrates that 
the asymptotic becomes discontinuous at b = 2Ra in the limit of large A. 



5 Wrong ways to study the asymptotic 

In this section we illustrate how unwarranted applications of the saddle point method 
can easily lead to incorrect results for the asymptotic of the tree approxumation to the 
Glauber AB amplitude, in particular to the results found in refs. [TJ and [2]. 



One may try to study the asymptotic directly from the representation (112p without 
passing to variables u and v. For A = B it can be rewritten as 

Q(h) - A ^ I dTdT ' J(1-Wt+t') A±(uv+u+v) ( o A \ 

where w and w are determined via r and r ; by Eqs. (11 II) and k and A are defined by 
(TTUT) and (1291) respectively. The results of pQ are obtained if we separate a factor in the 
integrand 

e «(T + T')-Aln(TT') 

and consider it as the only rapidly changing one at A — > oo. Then the saddle points are 

IT = it' = A. 

Putting these values in the integrand one obtains the asymptotic of the S'-matrix in the 
form 

S(b) = e F ^ b \ 

where 

F(b) = —iuv + u + v)- 2 AX 

K 

and u and v are determined by 

u = —iAne~ v , v = —iAne~ u . 

This is the result of pQ for the case of A = B and constant profile functions. 

However this derivation is too crude, since it neglects the A dependence of the ac- 
tion. A more elaborate derivation based on variables u and v leads to a pure optical 
approximation. Introducing for A = B integration variables t and t' defined as 

it = At, it' = Bt' 

we find 
where 

Pi = - ln(tt') + (t + t') + -(vu + u + v - r - r') 

7 

with r = yt and r' = jt', u and v determined via t and t' by equations 

u = rV, v = re~ u (36) 

and 7 defined by Eq. (j!7j) . Applying the saddle point method to the integral fl35|) one 
searches for saddle points from equations 

dPi = m = 

dt dt> 

Elementary calculations using Eqs. (IHUj) give 

dP x _ l \v dPx _ 1 Am 



Id 

If additionally A = 1 (central collisions) then the first and third terms cancel and the 
saddle points are found to be 

u = v = 7. 

They are the same as we had earlier, after passing to variables u and v. They lead to the 
optical approximation (see (1251) ). However we do not find the other pair of saddle points 
u — v — — 1. The reason is that in variables t and t' this saddle point is transformed into 
a singularity in the t, t'-plane present in the solutions of Eqs. flHol) . This singularity takes 
the leading role in the asymptotics at 7 > 1 and leads to the growth of the 5*-matrix and 
violation of unitarity. 

6 Conclusions 

Using the simplified form of the profile functions, constant inside the colliding nuclei, we 
have found that the set of tree diagrams for the AB amplitude in the Glauber approxima- 
tion is not unitary for heavy particpants and realistic values of the NN cross-section. This 
fact has been found analytically, using the saddle point method in the adequate manner, 
and fully confirmed by the straightforward numerical calculations. Unitarity is found to 
be fulfilled only if the NN cross-section a is small and diminishes with A (see Eq. (|2(jp ). 
Previous optimistic results [H [2] are found to be incorrect due to inadequate application 
of the saddle point method. 

From the practical point of view our results mean that the treatment of the AB 
scattering in the standard Glauber approximation must inevitably include loop diagrams 
and so depend on the form of the differential NN cross-section at non-zero angles. It 
remains to be studied how the situation changes if AB scattering is not described by the 
Glauber formula but is rather formed by the exchange of self-interacting pomerons (like 
in the Regge-Gribov model). In this case loops can also be formed, which are usually 
neglected because formally they are subdominant in the parameter A^ 1 ^ 3 . However it is 
not clear if the resulting tree amplitude is unitary when the limitation on the number of 
interacting nucleons at finite A is correctly imposed. We leave this problem for future 
studies. 
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8 Appendix. Derivation of the integral representa- 
tion for S-matrix. 

Consider the contribution to the AB scattering amplitude with a given number c of 
connected parts and a given number of participant nucleons n.j (n£) from the projectile 
(target) in the i-th connected part (i = l,...,c). The standard derivation leads to the 
expression 

iA c , niK (b) = iV—— — - -(«,)' J] / tPxaZMTpixi - b). (37) 

ci(/i — ny\n — n y. i=1 J 



1-i 



Here / is the total number of interactions: I = J21=i h where U is the number of interactions 
in the z-th connected part. Likewise n = Y%=i n i an d n' = Yh=i n [ are the total numbers of 
participants in the projectile and target. In the tree approximation we study Zj = rij+n^ — 
1, so that the number of interactions is uniquely determined by the number of participants. 
The symmetry factor N arises because with a given number of participant nucleons in 
each connected part there may be several terms which give identical contribution. To find 
this factor one may consider an auxiliary zero-dimensional quantum field theory with a 
generating functional 

Z = J D<pD<fte iY ^'^ ,p,p '\ 

where 

Y = a'VV - ip(e^ - l) - ip/U^ - l). 

This theory will generate the same diagrams as the Glauber expression fl37|) except that 
coordinate dependent densities Ta(x) and T B (x' — b) will be substituted by constants 
p and p' . The sum of all connected diagrams will be given by the effective action W 
obtained after integrating out the fields <fi and <jy. It will depend on powers of p and p' 
corresponding to numbers of fields entering in different connected parts 

i 

The coefficients C nji „' are just the symmetry factors which should be taken into account 
in Eq. ( 1371) except that they also include the corresponding number of amplitudes ia. So 
we obtain 

iA cmK {b) = d ( A _ n y( B _ n ry U C ni,n' i J ^x^^T^ix'*). 

Since we are interested only in tree diagrams, the effective action W should be taken in 
the classical approximation: 

W(p,p>) = Y( ( j ) (p,p'), ( J ) \p,p'),p,p'), 

where the classical fields are determined from a pair of transcendental equations 

(f) = iap / e </ '\ <f) !f =iape^. 

Next we sum over all posible values of n«, n- and c. To sum over rij and n! i we present 

Summations over rij and n! i in the integrand factorize in two factors 

A 1 c c 
m l A _ n >- i=i 1 

where Ui = Tj^ix^jzi and a similar factor for the target B. We present each as 

1 f°° i 
u n = — dtt n e- t,u 
un\ Jo 



1-J 



to find a sum in the integrand 

Y / a ' \i TT = ( 1 + Y U) - zero terms. (38) 

The subtraction terms eliminate contributions from the z-th connected part when all 
rii = 0. The right-hand side of Eq. (1381) can be represented as a contour integral 



A\ r dr 



i=i 

Integration over all ij gives 



. , A\ f dr , 1 

/A( ^ ) = 2^/^niT3— -- 1 



The final expression for the amplitude is obtained after integrations over all Zi and sum- 
mation over the number of connected parts c, which are realized in a straightforward 
manner. Changing <fi — > —u and <f>* — > — u leads to the expression Eq. (jSJ). 
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